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Abstract. 

We deal with a model for high-temperature superconductivity which main- 
tains that in cuprates electrons running in the copper oxide layers, found in 
lattice of these materials, form spin-singlet bonds with electrons running in the 
neighbouring layers. This model reutilizes the BCS scheme, but with the es- 
sential difference that the electron pairs are characterized by equal, rather than 
opposite momenta as in Cooper pairs. In the present paper, we consider the 
electron pair formation and a peculiar canonical transformation analogous to 
the transformation once applied to the theory of pairing correlations in nuclear 
matter. It is shown that the quasi-particle energy spectrum remains that of 
the BCS theory, including the linear relationship between forbidden energy gap 
and critical temperature. The model is also applied to superconductivity of 
some copperless perovskites of mixed stoichiometry, whose features are of spe- 
cial worth in understanding the mechanism of the phenomenon. The possibility 
of enhancing critical temperature in cuprates by inserting monovalent ions into 
the lattice is considered. 

PACS: 74.20.-z; 74.72. JT. 
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1. Introduction 

Notwithstanding the great deal of work done till now, no theory about high 
temperature superconductivity has obtained a general consensus. It is common 
belief that even the basic nature of the phenomenon is not understood. In maga- 
zine notes appeared on "Scientific American" in 2000 and 2004, the inadequacy 
of theoretical models in explaining superconduction in cuprates is stressed [1] . 
Owing to this state of affairs, we will now examine a mechanism, quite unlike 
those so far proposed, which has been conceived by keeping in mind, besides 
cuprates, the features of other kinds of unconventional superconductors different 
from cuprates ( 2 ). 

Cuprates are surely the most interesting superconductors as they allow for 
the highest critical temperatures so far recorded. But a variety of materials 
other than cuprates is known, showing a superconductivity not explained by 
the BCS theory. Actually, superconduction has been detected in perovskites 
of fractional stoichiometry, in mixed copper and alkaly-earth oxides, in organic 

1 Corresponding author E-mail address: brovetto@vaxcal.unica.it. 
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2 ) In Ref. [2, 3, 4, 5] , some features of this mechanism have already been presented. 
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compounds and in fullerenes. This makes complex the study of superconduc- 
tivity but, at the same time, provides a wide experimental basis with which 
theory must be compared. In our opinion, the simplest and most conservative 
hypothesis is that the essential superconduction mechanism is the same in all 
these materials in spite of their quite different natures. Accordingly, the very 
cause of superconductivity must be searched for in something which is shared 
by all materials. On that account, next Section is devoted to singling out the 
features which pertain to all above-mentioned superconductors. 

2. Common features of unconventional superconduc- 
tors 

In reality, two remarkable peculiarities are shared by the superconductors 
cited before. The first is that all are characterized by complex layered lattices 
or uneven heterogeneous lattices showing discontinuous structures. The second 
is that all contain ions or atoms with unpaired electrons or electrons not in- 
cluded in closed shells. These points are emphasized hereinafter by considering 
a selection of various superconductors. 

A ) Perovskites with fractional stoichiometrics. - Some of these perovskitcs 
are listed in the following Table. 

Table 1: Superconducting perovskites with fractional stoichiometries. 

SrTi0 3 -^ [6] BaPbo.7Bio.3O3 [7] Ba .6Ko.4Bi0 3 [8] 
T C ~0.3K T C = 13K T C = 30K 

Sr .5Ko.5Bi0 3 [9] Sr . 5 Rbo.5Bi0 3 [9] 
T C = 12K T C = 13K 

Their structure is characterized by lattice discontinuities found at the bor- 
ders between cells with different ion compositions. The SrTi03_,5 supercon- 
ductor, which shows a partial lack of oxygen, is a reduced compound. But, 
also BaPbo.7Bio.3O3 and Bao.6Ko.4Bi0 3 are indeed reduced compounds. In 
fact, owing to valence four of lead and five of the fully oxidized bismuth, their 
stoichiometries should be written as BaPbo.7Bio.3O3. 15 and Bao.6Ko.4 Bi03. 3 , re- 
spectively. The same argument obviously is right for the strontium-substituted 
compounds. For these compounds, the lack of room in the stiff perovskitic cell 
prevents oxygen from entering the cell until metals are fully oxidized. Since oxy- 
gen is kept in the form of divalent 0~ 2 ions, when oxygen is removed as neutral 
atoms some electrons are left in the material and become bound to metal ions. 
It follows that unpaired electrons appear in excess to the noble gas shells of 
K+ 1 , Ba+ 2 , Sr+ 2 ions or to 5d 10 shell of Pb+ 4 and Bi+ 5 ions. 

B) Cupratcs. - These materials show layered lattices formed by perovskitic 
or perovskitic-like cubes. As cuprates are the best known superconductors, we 
limit ourselves to few examples. The first discovered Lai.85Sro.i5Cu04 cuprate, 
which superconducts at 35 K, is characterized by the K2MF4 structure, that is, 
an alternation of perovskitic and NaCl-like layers [10] . It shows a fractionary 
stoichiomctry and is to be regarded as an oxidized superconductor because, 
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owing to substitution of trivalent lanthanum with divalent strontium, its stoi- 
chiometry should be written as La1.g5Sro.15CuO3.925. The 92 K superconductor 
YBa 2 Cu 3 07, usually referred to as YBCO, is characterized by a stacking of yt- 
trium and barium centred lacunar perovskitic cubes [11] . In the so-called TBCO 
superconductors, such as for instance the Tl2Ba2CaCu20s compound, different 
alternation of perovskitic-like layers of copper, calcium, barium and thallium 
oxide are found [12, 13] . As for the unpaired electrons, we point out that all 
cuprates contain divalent copper with the [Ar]3d 9 configuration showing just 
one unpaired 3d-electron. 

C) Mixed copper and alkaly-earth oxides. These compounds deserve atten- 
tion because they are cuprates lacking in perovskitic structure. The mixed 
oxide SrCu0 2 shows an orthorhombic lattice, it is not a superconductor but 
superconductivity appears at 40 K in the fractionary stoichiometry compound 
Sro.86Ndo.i4Cu02 [14] . It is a reduced compounds because, owing to valence 
three of neodimium, its stoichiometry should be written as Sro.86Ndo.14CuO2.07- 
Apart from the different structure, it is like the superconductors of item A ). Re- 
cently, using a field-effect technique, electrons were removed from (or injected 
into) the monoclinic CaCu02 compound. In this way, superconductivity was 
found at 89 K and 34 K depending on wheter 0.15 electrons per molecule are 
removed or injected, respectively [15]. Even in this case, superconductivity is 
originated by introduction of unpaired electrons and of lattice discontinuities 
lying at the borders between cells of different degrees of oxidation ( 3 ). 

D) Organic superconductors. We limit ourselves to the Bechgaard salt, that 
is, tetramcthyl-tctraselena-fulvalcnc hcxafluoro-phosphate (TMTSF)2 PFg which 
superconducts at about 1 K [16] . This material is characterized by stackings of 
strongly bound molecules with a much weaker intermolecular bonding in di- 
rection transverse to the stackings. One electron is moved from one TMTSF 
molecule to one fluorine atom so that (TMTSF) +1 cations and (PFg) -1 anions 
appear. Since in the neutral TMTSF molecules all electrons are coupled in o— 
or tt— bonds, one unpaired electron is present in the (TMTSF) +1 cation. 

E) Fullerenes. These materials are characterized by stacks of Cgo balls. Links 
between carbons in contiguous balls are weaker than those of carbons in the 
same ball. Superconductivity has been detected at 18 K in K3C60 and at 28 
K in RbaCeo [17] . The presence of potassium or rubidium atoms, showing the 
[Ar]4s and [Kr]5s configurations, respectively, inserts unpaired electrons in the 
Ceo stacks. With the field-effect technique, electrons were removed or injected 
into the Ceo balls, so leaving some unpaired electrons there. In this way, super- 
conductivity was originated at peak temperatures of 52 K or 11 K when just 
three electrons were taken off or added to each ball, respectively [18] . 

The previous analysis confirms that uneven lattices and unpaired electrons 
are really features common to the superconductors considered above. How- 
ever, different kinds of superconductors must be distinguished depending on 
the actual provenance of the unpaired electrons. Indeed, compounds of item 

3 ) With the field-effect technique the average degree of oxidation of the material can be 
properly determined. On the contrary, lattice discontinuities related to the local degree of 
oxidation of the cells remain uncertain. Also incidental lattice defects might play a role. 
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A) and the Sro.86Ndo.i4Cu02 compound are "reduced" superconductors. The 
Lai.85Sro.i5Cu04 cuprate, on the contrary, is an "oxidized" superconductor. 
Cuprates as YBa 2 Cu 3 07 or Tl 2 Ba 2 CaCu208 and the Bechgaard salt are to be 
regarded as "intrinsic" superconductors, since unpaired electrons are peculiar 
to their chemical composition. The alkaly-doped fullerene as well as fullerene 
and CaCuC>2 oxide showing field-effect superconductivity are to be regarded 
as "doped" superconductors, because in these materials superconductivity is 
originated by an external agent. 

In our opinion, the features examined above can be considered as a sort of 
"Rosetta stone" for disentangling the problem of high temperature supercon- 
ductivity. 

3. About properties of fermion systems 

Let us recall some topics concerning properties of fermion systems which 
will be helpful in understanding the mechanism of superconduction in the ma- 
terials cited above. In 1916, G.N. Lewis first discovered that covalent bonds 
consist of pairs of shared electrons [19]. This fact, inexplicable by the classical 
physics, was interpreted in 1927 by W. Heitler and F. London (HL) who ap- 
plied quantum mechanics to the hydrogen molecule [20]. By considering two 
hydrogen atoms A and B in Is states, they wrote a two-electron wave func- 
tion of the form: [ui s a (1) u\ s b (2) + ui s b (1) mi s a (2)] in which each electron 
is found at the same time both on atom A and B. This function, symmet- 
ric with respect to exchange of electrons, was associated to an antisymmetric 
spin function: [a (1) (3 (2) — (3 (1) a (2)] representing a spin-singlet state, so al- 
lowing for the Pauli principle. In this way, in evaluating the expectation value 
of energy, integrals involving products of electron states: u\ s a (1) uub (1) and 
u\sA (2) u\sb (2) appear in calculations. These exchange integrals account for 
covalent bond energy. A year later, W. Heisemberg, utilizing the same argu- 
ments, explained the origin of the Weiss field in ferromagnetic solids [21]. 

In 1933, exchange forces came back into evidence in a quite different field of 
physics. In this year, indeed, E. Majorana, in dealing with nuclear interactions, 
introduced forces which exchange the coordinates of the interacting nucleons 
[22]. These forces are mediated by charged pions and act only for nucleons in 
neighbouring momentum states [23]. The Majorana forces are of paramount 
importance since they account for the saturation effects in binding energy of 
nuclei. 

In 1957 the famous BCS theory finally explained superconduction in metals 
[24]. The essential device of this theory are the Cooper pairs, that is, pairs of 
electrons of opposite momenta bound by a phonon coupling. Utilizing a special 
canonical transformation devised by N.N. Bogolyubov, the system of interacting 
electrons is substituted by a set of non-interacting quasi-particles showing an 
energy gap at the top of the distribution [25]. 

The great success of the BCS theory drew attention on the possibility of 
its application to nuclear physics. In 1958, A. Bohr, B.R. Mottelson and D. 
Pines proposed that the energy gap found in the spectra of even-even nuclei is 
originated by a mechanism analogous to that of superconduction in metals [26] . 
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Indeed, when the main part of the nucleon-nucleon interactions is averaged so 
allowing for a self-consistent field, there remains a residual weak interaction, 
due to Majorana forces, which couples nucleons of like momenta. A thorough 
treatment of this problem was performed by S.T. Belyaev which modified the 
Bogolyubov transformation substituting the Cooper pairs with Majorana pairs 
of nucleons with equal linear momenta, but opposite projections of angular 
momenta along the quantization axis [27] . This treatment, however, leaves out 
the dependence on temperature of the energy gap, owing to the fact that nuclei 
are always on the ground state. An equivalent treatment was performed by L.P. 
Gor'kov and A.I. Alcksccv utilizing the Green function technique [28, 29]. 

Since our treatment on superconductivity utilizes a method similar to that 
applied by Belyaev, an extensive account on this matter is given in Appendix 
A. 

4. The superconducting Lewis pairs 

The superconductor features highlighted in Section 2 and the arguments 
presented in Section 3 induce us to argue that at low temperature unpaired 
electrons running in a superconductor region bordering on a lattice discontinu- 
ity originate spin-singlet pairs with electrons running in the region bordering 
on the opposite side of the lattice discontinuity. This is due to instability of the 
unpaired electrons that tend to form covalent bonds. Obviously, in order to set 
out a quantitative treatment, it is necessary to know, besides the electron wave 
functions, the actual nature of the lattice discontinuities. This occurs with the 
intrinsic superconductors, such as the the YBa2Cu307 and the Tl2Ba2CaCu20s 
cuprates, or the Bechgaard salt. On the contrary, with the reduced supercon- 
ductors of items A ) and C) it is necessary to resort to special conjectures, since 
their fractional stoichiometries make the structure uncertain. These supercon- 
ductors are investigated in Section 7. Also the fullcrcnc-based supercondutors 
give rise to difficulties of this kind. In practice, the materials most right for our 
investigations are the intrinsic cuprates. 

The previously cited cuprates, are characterized by planes of oxygen lacunae 
and yttrium or calcium ions sandwiched between couples of contiguous Cu0 2 
layers (see [30] Ch. 7). In the following, these layers will be marked with labels 
a and b. Two unpaired electrons, one running on layer a the other on layer 6, 
can be represented by the tight-binding (TB) wave functions 



in which k a and kt mean the electron wave vectors, a(ri — u p ) and 6(r 2 — 
v q ) the 3d-orbitals of the copper ions on layers a and b and u p and v q their 
lattice vectors, respectively. Each copper ion on layer a is separated from a 





(1) 
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corresponding ion on layer b by the spacing A between the layers, that is, 

v p - u p = A. (2) 
The energies of the unpaired electrons spoken of are 

W a (k a ) = (0 Q (k Q ,n)| H a (pi,n) |0 a (k Q ,n)) , 

W b (k b ) = (fa (k b , r 2 ) | H b (pa , r 2 ) | & (k 6 , r 2 )) , (3) 

where the Hamiltonians -ff a and Hb account for the electron kinetic energies 
Pi/2m and p 2 /2m and for the Coulomb interactions of electrons 1 and 2 with the 
copper ions in positions u p and v q , respectively. For k Q = k b , energies W a (k Q ) 
and Wb (kb) are equal, owing to equality of the Cu0 2 layers. In Appendix B, 
utilizing a special model for the actual nature of the copper ion orbitals, energies 
W a (k a ) and Wb (k b ) are evaluated on the ground of Eqs. ©. 

Owing to the peculiar structure of the before cited cuprates, that is, the 
presence of oxygen lacunae on the yttrium or calcium planes placed between 
the copper ions, an unpaired electron of layer a is allowed to form a covalent 
bond with an unpaired electron of layer b, like the Is electron of a hydrogen 
atom A forms a covalent bond with the Is electron of another hydrogen atom 
B. On this ground, in analogy to the HL treatment of the hydrogen molecule 
[31], the wave function for a pair of electrons of layers a and & in a spin-singlet 
state is 

* <>i, r 2 ) = 1 ===== [<f> a (K, ri)</> b (k b , r 2 ) + Q (k Q , r 2 )0 b (k b , n)] x 

x-L[a(l)/3(2)-a(2)/3(l)], (4) 

a (1) and (3 (2) standing for the spin functions. In the following, these pairs 
are referred to as "Lewis pairs" since this author, already cited in Section 3, 
pioneered investigations on covalent bonds. The possibility of applying the HL 
treatment is due to the fact that it is implemented aside from the actual nature 
of the electron states, so that ls^ and Isb or 4> a and 4> b states can be indifferently 
considered. This notwithstanding the fact that Is states account for a single 
Coulomb potential centre, while 4> states account for N centres. Like in the HL 
treatment, energy —Wp of the electron pair holds a "classic" contribution, that 
is, without exchange of electrons between (j> a and </> b states and an exchange 
contribution in which both electrons are shared between 4> a and <fi b states ( 4 ), 
that is, 

-Wp = (* (ri ,r 2 )|g int (ri,r 2 )|tf (r^ra)) = 

4 ) Apart from substitution of Isa and lsg hydrogen-like states with <f> a and <j> b states and 
the presence of summations over the copper ions, the terms appearing in Eq. (5) are like 
the terms in Eq. (43-7) and (43-9) of reference [31] p. 342. We omit considerig terms for the 
spin-triplet state which, in HL treatment, originate repulsion between the hydrogen atoms. 
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In this equation, Hamiltonian i?i„t (ri,r 2 ) allows for Coulomb interactions 
of electrons with copper ions of effective charge Z and for Coulomb repulsion 
between the electrons. Like in the HL treatment, it follows that pairing energy 
is given by 
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(6) 



Terms J and J' represent the classic contributions to energy of the electron 
pair. These contributions are negligible as occurs in the case of the hydrogen 
molecule [31]. The squared overlap integral (<fi a | b ) 2 is negligible with respect 
to unity. Only integrals K and K', which account for exchange interactions, 
must be retained. 

We proceed now to evaluate the integrals appearing in Eq. |@. Let us first 
consider the exchange integral K' . Taking into account Eqs. QJ, we have 



K' = <0 a (k a , n) b (k b , r 2 )| \<j> b (k b , n) a (k a , r 2 )) 
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d 3 r 2 . (7) 



This equation involves sums over TV 4 terms. But, taking into account that the 
electron distributions in orbitals a(ri — u p ) or a(r 2 — u p ) and fe(ri — v g ) or 
6(r 2 — v g ) are closely localized at the lattice positions u p and v g , respectively, 
terms b (r 2 — v g ) a (r 2 — u p ) and a (r\ — u p ) b (ri — v q ) are non-negligible only 
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when T2 — Vg and r 2 ~ u p and ri ~ u p and ri ~ v g , which entails that v g ~ u p , 
that is, remembering Eq. (J2J, p = g. It follows that the sums in Eq. J7J) contain 
only N 2 non-negligible terms. This leads to 



e 2 f A 

K ' = z2 exp I* ( ka ~ kf> ) K ~ u p ~ A )i x 

g,P=i 



x 6 (r 2 - a (r 2 - u q ) a (n - u p ) 6 (n - v p ) (8) 

r l,2 

In this sum, significant contributions appear only when n,2 is small, that is 
when r2 — ri ( 5 ). This, as before, entails that only terms with v q ~ u p , that 
is, p = g are non-negligible. So we obtain 

K' = — f V & ( r 2 ~ V <?) a ( r 2 ~ U <?) Q ( r l ~ U g) b ( r l - V g) .3 

iv 2 y £j n, 2 1 2 



(9) 

which means that this integral, which accounts for Coulomb repulsion between 
electrons, can be neglected. We consider now the overlap integral 

(cj) a (k„,ri) | 4> b (k b ,ri)) = 



I f N 

= ^2 exp(-ik a -u p + ik b -v q )a(r 1 -u p )b(r 1 -v q )d 3 r 1 . (10) 

J p,q=l 

In this case, significant contributions are obtained only for orbitals with u p ~ ri 
and v q ~ ri which entails that u p ~ v q and, as before, p = q. We have thus 

(<j) a (k a ,ri) | (j) b (k 6 ,ri)) = 

f 1 N 

= / exp [-i (k a - k 6 ) ■ r Ml ] — a (n - u p ) 6 (n — u p — A) d 3 ri, (11) 

in which only the parallel component r Ml has been accounted for in the expo- 
nential factor since k a and kf, are parallel to a and b layers. By putting 

i N r 

F ^ = wE / a (r x - Up) 6 (rj — u p — A) rfr ± u (12) 



5 ) It is to be pointed out that Coulomb repulsion between electrons is screened by the 
presence of the positive ions which assure the electric neutrality of the material. Consequently, 
Coulomb repulsion decreases more quickly than e 2 /ri 2- This is like what occurs with Cooper 
pairs in BCS theory. 
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Eq. qilfl can be rewritten as a Fourier transform, that is, 

{<f> a (k OJ n) I (f> b (k b , t x )) = J exp [-« (K - k b ) ■ r„i] F (r llX ) d 2 r nl . (13) 

Function F (r„i) shows a slight periodic dependence on r„i since, when r„i varies, 
about equivalent terms are always included in the sum over p. In particular, if 
dependence on r M i of function F (r M i) is omitted we have 

<0 Q (k Q , ri ) | <S> b (k b , n)) = (2^) 2 F <5 (k„ - k b ) . (14) 

Eqs. (|13|l and (|14|l entail that electrons of equal wave vectors, that is, of equal 
momenta, allow for the maximum value of overlap integral and, therefore, the 
maximum pairing energy. It is worth to point out that this result mimics that 
concerning the Majorana exchange forces, which likewise originate fermion pairs 
of equal momenta. But, with the Lewis pairs this is a consequence of antisym- 
metry of the wave function J3J which accounts for the Pauli's principle, while 
with the Majorana pairs equality of momenta directly follows from the exchange 
nature of the forces, which are mediated by charged pions. On this ground, we 
consider in the following only electron pairs with k a = k b . We find in this way 



1 N f 

{<t> a (k.,ri) | cj> b (ki,ri)) =mY1 / a ( r i- u p) & ( r i- u p--V>d 3r i = S a< b, (15) 

P =i J 

where 

S a , b = J a(p)b(p-\)d 3 p. (16) 
By applying the same procedure to exchange integral 



K = - (^(k 6 ,r 2 )| T^—\ I6,(ka. r 2 )) , (17) 
P =i l r2 - u fl 



we obtain 



I N r 2 2 
K =~TfYl / a ( r 2- u 0i = j-6(ra-uj- \)d 3 r 2 =E b , a , (18) 



l=i 

where 



f Ze 2 

E b>a = - / a(p)- -b(p-\)d 3 p. 

J \P\ 



(19) 



Consequently, utilizing Eqs. ©, ©: <|15[) and l|18fl . pairing energy for k Q = k b 
turns out to be 

-W P = 2 5 M £ M . (20) 
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Energy Wp, of course, is expected to be very small in comparison with those 
of most covalent bonds, owing to the considerable separation of the interacting 
copper oxide layers which is larger than normal covalent bond lengths. 

5. The canonical transformation 

According to the above picture, electrons running on the contiguous copper 
oxide layers a and b constitute one set of 2N Fermi's particles. An electron in this 
set with a given wave- vector k shows degeneration due to its spin components a 
and P and degeneration due to options k = k a or k = k b , that is, to its placing 
on layer a or b. Owing to k degeneration, the electron set is split into two 
conjugated subsets a and b. Each electron of subset a is paired with an electron 
of subset b of equal wave vector, that is, of equal momentum, with an energy 
—Wp independent of the actual momentum. As follows from spin-singlet wave 
function (J2J| , these pairs show spin zero so that can be regarded as weakly bound 
bosons. On this ground, the second quantization Hamiltonian for Lewis pairs 
can be identified with the Belyaev's Hamiltonian 14511 considered in Appendix A 
in connection with systems of interacting particles characterized by symmetry 
of reflection in a plane. This, of course, barring substitution of indexes R and L 
with indexes a and b. In the case of YE^CuaOy cuprate, the plane of symmetry 
is to be identified with the plane of yttrium ions. In this way, by taking single 
particle energies into account, writing briefly s k for W a (k Q ) and W b (k & ) and 
substituting Vk,k> with Wp, the Hamiltonian for a system of electrons forming 
Lewis pairs can be written in the quick form 

H = ( £fe ~ V) ("fea"fca + a% b a kb ) -W P ^2 ut>a°i, b u kb a ka (21) 

/j, standing for the chemical potential of the 2N electron set. In Appendix 
C, simple arguments are posed showing that this Hamiltonian conserves the 
momentum of the system. 

Assuming coefficients which satisfy condition 

U 2 k +V k 2 = l, (22) 

the canonical transformation is 

"fca = U k (3 ka + V k kb , 
a kb = U k p kb - V k /3 ka , 

P ka and [3 kb standing for the quasi-particle destruction operators. Substitution 
of Eq. (23 into Eq. J2TJ leads to 

H = J2 Z k [%V k 2 + (t/f - V k 2 ) (n ka + n kb ) + 2U k V k ffijtb + KbPka) ■ 
k 

~WpJ2 S A ( 24 ) 

fc, k' 
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where £ fc = e k — \i and 

Bk = Ufp kb p ka - V k 2 p + k X b + U k V k (1 - n ka - n kb ) , (25) 

n ka = (3 ka (3 ka and n kb = f3 kb /3 kb standing for the operators of quasi-particle 
numbers. These equations, apart from some obvious differences in symbols, are 
like those for the BCS theory given, for instance, in the Landau and Lifshitz 
treatise [32]. For this reason, reutilizing the same routine procedure, we limit 
ourselves to reporting the most significant issues. By taking condition ill' I'll into 
account, considering the diagonal terms in the Hamiltonian and minimizing 
energy with respect to coefficient U k for given entropy, we find 

2£ k U h V k = A (U 2 k - Vi) , (26) 

with 

A = W P J2 U k'Vk'(l-n k > a -n klb ), (27) 
fe' 

in which Ti k i a and n k >b now mean the actual numbers of quasi-particles. When 

Eq. 126fl is verified, the non-diagonal second-order terms P kb P ka and j3 ka (3 kb are 
removed from Hamiltonian l)24[l. From Eqs. 1(22(1 and 1)26(1 . the usual relations 
for coefficients U k and V k are obtained, 
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which, when substituted into Eq. 1(27(1 . lead to the equation which determines 
A 



The next step is to change summation on k' to integration on energy. Taking 
into account that for e k — and e k = we have £ fc = — /U and £ fe = 0, 
respectively, Eq. ((29(1 becomes 

w P f (1 ~^ a ~"f ) « (o = i, (so) 

fi(£) standing for the density of states per unit cell. Owing to the expected 
smallness of A with respect to /i, the main contribution to the integrand arises 
for £ ~ so that we can put Q (£) ~ Qp. For £ = 0, we have indeed £& = fi = £p 
since the chemical potential, apart from a small correction due to temperature, 



11 



coincides with the electron kinetic energy at the Fermi level ( 6 ). Consequently, 
at T = K where n ia = n^ b = 0, Eq. (JSDJl yields 



which leads to ( 7 ) 



WpVlF log 



A 



v/Al+7^ 



W>^log^ = l (31) 



A = 2/i cxp (-i/w P n F ) ■ 



(32) 



This equation differs from the corresponding one of the BCS theory only in the 
substitution of the Debye energy Hlud with twice the chemical potential fi. It 
follows from Eq. (|32|l that superconduction is ruled by three parameters, that is, 
pair energy Wp, density of states Oj? at the Fermi level and chemical potential 
fi. Taking into account that /i = sp, a connexion between /i and Q F is expected, 
depending on the actual electron energy spectrum ( 8 ). In Appendix B, utilizing 
special assumptions for the copper ions orbitals, this connexion is found to be: 
fiflp = 1.28, which allows Eq. Ilii2[l to be rewritten as 



(-0.78^ 



For T > 0, substituting 



n (a = n ib = 1/ 



exp 



A 2 + £ 2 /kT + 1 



into Eq. JSUJl, we obtain as in the BCS theory 



log 



7Ao 7C(3) 



A 



nkpT 8tt 2 \ksT 



(33) 



(34) 



(35) 



where kp is the Boltzmann constant, log 7 = 0.577 the Euler constant and 
( (3) = 1.202 the Riemann Zeta-function. By putting A = in this equation, 
the usual relationship between energy gap and critical temperature is found 



2A = — k B T c = 3.52k B T c 
7 



(36) 



6 ) See for instance [33] Ch. III. 

7 ) In this equation Ao is a constant quantity. But, if in normalizing the TB functions iTi . 
overlap integrals of copper ion orbitals are not disregarded, Ao is substituted by a quantity 
A (9) depending on the angle between the electron wave vector k and the cell a-axis. This ac- 
counts for the d-symmetry of the order parameter. In Ref. [3] this matter has been thoroughly 
discussed. 

8 ) For a tridimensional Fermi gas of TV electrons in a volume V, we have: NQp = 



(V/27T 2 ) (2m/h 2 ) 



3/2 



where: 



leads to: fiflp = 1.5 (see for instance [33] Ch. III). 



(h 2 /2m) (3n 2 N/V) 



2/3 



This, by letting /i = ep, 
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Likewise, by utilizing the Hamiltonian (|24[1 and Eqs. (|27l) and l|28|) . the quasi- 
particle energy spectrum is found to be the same as that of the BCS theory, 
that is, 

w k = ^A 2 + (e fe -e F ) 2 . (37) 

The actual magnitude of the energy gap is 2Ao since quasi-particles appear in 
pairs as occurs in BCS theory. It is to pointed out that Eq. i|36[) has been 
successfully tested for various unconventional superconductors. An extensive 
tabulation of data concerned is given in [30] Ch. 6. 

6. Experimental evidence in favour of the interacting- 
layer mechanism 

We will now briefly examine some experimental results which substantiate 
the interacting-layer mechanism. They are reported here in order of increasing 
significance. 

1 ) Coherence length - A first clue about the interaction between the super- 
conducting layers is offered by measurements of the Hall effect on YBa2Cu307 
single crystals. It was found that the Ginzburg-Landau coherence length along 

o 

the c-axis is 1.5 A. Since the spacing of superconducting layers is near to c- 

o 

axis lattice parameter 11.68 A, the conclusion was drawn that "t/ie two copper 

o 

oxide planes which are spaced 3.2 A, are tightly coupled and act as a single 
superconducting layer" [34] . 

2) The effect of internal pressure. - More direct evidence comes out from 
the effect of pressure on critical temperature. It is common knowledge that in 
cuprates T c considerably increases when samples arc submitted to hydrostatic 
pressures on the order of few GPa. Eqs. (|16|) and (|19|) . which relate pairing 
energy to the layer separation A, explain this effect. Indeed, hydrostatic pressure 
lessens separation A thus increasing Wp and, consequently, critical temperature. 
But hydrostatic pressure lessens the distances between copper ions in direction 
both parallel and orthogonal to CuC>2 layers. The effect of hydrostatic pressure, 
therefore, is unsuitable in distinguishing interactions inside each layer from those 
between contiguous layers so that no evidence in favour of the interacting-layer 
mechanism is obtained. The conclusion, however, is different if the so-called 
"internal" or "chemical" pressure is considered. Indeed, substitution of some 
ions with others of smaller radius originates a decrease in the cell size which is 
commonly regarded as the effect of an internal pressure. In this connection, let 
us quote, the following sentence by P. Chu highlighted in a note by K.A. Muller 
[35] : " Therefore, Paul Chu thought, O.K., instead of applying pressure, I rather 
use a rare-earth ion, namely the yttrium which is smaller than lanthanum, and 
thus get a higher T c owing to the induced internal pressure " . Actually, the La + 

o o 

ionic radius is 1.15 A while that of Y +3 is 0.93 A. This is tantamount to saying 
that in YBCO T c increases just when the contiguous CuC-2 layers approach each 
other leaving unchanged copper ion distances parallel to the layers ( 9 ). 

a ) A clear evidence of the effect of internal pressure is offered by some thallium-based 
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3) The effects of yttrium-praseodymium substitution and of reduction in 
YBa 2 Cu 3 07. - Concerning the effect of ion substitutions in YBa 2 Cu307, a sur- 
prising feature is the lack of superconductivity of the PrBa 2 Cu 3 C>7 compound 
[37, 38] . This fact cannot be ascribed to the praseodymium ionic radius which 
is not significantly different from the yttrium radius. To explain this peculiar re- 
sult it should be taken into account that praseodymium originates both trivalent 
and tetravalent ions. Actually, magnetic susceptibility measurements have indi- 
cated that, in the compound dealt with, praseodymium is tetravalent [39, 40] . 
This means that one electron is released from praseodymium and transferred to 
the neighbouring copper ions on the Cu0 2 layers so that each cell contains one 
divalent and one monovalent copper rather than two divalent coppers as in the 
yttrium compound. In the PrBa 2 Cu 3 07 compound praseodymium acts as an 
electron donor. Since monovalent copper shows the [Ar]3d 10 configuration lack- 
ing unpaired electrons, the superconducting pairs can no longer be originated. 
Consequently, superconductivity is shut out by the interacting-layer mechanism, 
just as expected. 

A very interesting result concerns the effect of reduction on YBa 2 Cu 3 07 
critical temperature. Samples of stoichiometry YBa 2 Cu307_a; show a decreasing 
T c for increasing x. Actually, the T c versus x plot is characterized by two 
plateaux, the first at 92 K, for x less than 0.2, followed by a step decrease 
and by a second plateau at about 60 K, for x near to 0.4. Measurements of 
distances of copper from nearby ions have shown that the effective valence ( 10 ) 
of copper on Cu0 2 layers is characterized by a parallel behaviour. Indeed, the 
plot of copper effective valence versus x shows just two plateaux for the same 
values of x separated by a step decrease of 0.03 e/Cu effective charges [44]. 
This behaviour, like that originated by the yttrium-praseodimium substitution, 
depends on reduction of divalent copper on the Cu0 2 layers. When oxygen is 
removed, a number of electrons is left into the lattice. Along the plateaux, only 
trivalent copper on the cell basal planes is reduced thus leaving T c unaffected. 
The T c decrease for 0.2 < x < 0.4 corresponds to the decrease of effective valence 
of Cu0 2 layer copper. In Reference [5] a thorough thermodinamic description 
of this effect is given. 

4) The monolayered Tl 2 Ba 2 Cu06 compound. - Several multilayered thallium- 
based superconductors are known. The compond mentioned here represents a 
special case since it is characterized by a single Cu0 2 layer interposed between 
two BaO and two TIO layers on the outside of the BaO layers (see [30] Ch. 7). 
With this structure, the Cu0 2 layers are well separated so that the interacting- 
layer mechanism cannot be active. In spite of this, this compound supercon- 
ducts at 85 K [13]. This fact may seem to represent strong evidence against the 
mechanism we consider. In reality, the situation is quite opposite. In fact, only 
reduced samples of stoichiometry Tl 2 Ba 2 Cu0 6 _5 superconduct. Experiments 

cuprates in which barium-strontium substitution slightly lessens the cell size along c-axis so 
originating T c enhancements as large as tens of K [36] . 

10 ) For the meaning of the "effective valence" parameter, otherwise referred to as the "bond 
valence sum", sec [41, 42, .43] . It can be identified, in practice, with the copper average 
valence. 
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have shown that in fully oxidized samples superconductivity is destroyed just as 
expected on the ground of the mechanism dealt with [45] . Reduction introduces 
in the lattice unpaired electrons, as occurs in the perovskites with fractional 
stoichiometry listed in item A) of Section 2. As it will be shown in the next 
Section, these electrons originate additional layers of unpaired electrons with 
which electrons in the CuC>2 layers interact so allowing for superconductivity. 
On this ground, the sudden onset of superconductivity observed as soon as the 
oxygen content is reduced is not surprising. Indeed, even a minimum quantity 
of superconducting material embedded in an inert matrix is sufficient to set the 
sample resistance to zero. 

By keeping the previous arguments in mind, the queer result that reduc- 
tion sometimes causes and other times hinders superconductivity is explained. 
In Ref. [4] , other evidences in favour of the interacting-layer mechanism are 
discussed. 

7. Superconductivity of mixed stoichiometry 
perovskites 

We extend our analysis to superconductivity of some mixed stoichiometry 
perovskites. As pointed out in Section 2, these materials are to be regarded as 
reduced compounds. Reduction decreases the actual cation valence, thus intro- 
ducing unpaired electrons into the lattice. In the mechanism we consider, un- 
paired electrons are indeed the basic ingredient for superconductivity. However, 
the question is to be settled of the lattice structure which allows the interacting 
layer mechanism to operate. While cuprates show a quite tidy layered structure, 
in the mixed stoichiometry compounds oxygen lacunae or substitutional ions are 
placed at random. Despite this, for mere statistical reasons it can be expected 
that the lattice contains some domains in which ions are layered in the right 
order to allow superconductivity. We point out, in this connection, that the 
formation of Lewis pairs is allowed even when lattice vectors v 9 and u p are not 
ordered on plane surfaces, as occurs in cuprates. Even irregularly bent surfaces, 
such as those that are likely found in mixed stoichiometry compounds, are suit- 
able, provided that a number of orbitals with unpaired electrons exist such that 
Eq. (J2J can be applied. Therefore, remembering that even a minimum quantity 
of superconducting phase sets sample resistance to zero, the superconductivity 
of mixed stoichiometry perovskites can reasonably be explained. 

To understand how this can occur, let us focus attention on the previously 
cited 30 K superconductor Bao.eKcuBiOa. The most conservative assumption 
is that the unpaired electrons introduced by the potassium-barium substitution 
lie right on the barium, so that monovalent Ba +1 ions substitute K +1 ions. 
In this way, in fact, the lattice Madelung energy is kept unchanged at its for- 
mer value. Evidence in favour of this assumption is offered by the fact that 
superconductivity was detected in the Sr .5K .5BiO3 and Sr .5Rbo.s BiC>3 com- 
pounds at 12 K and 13 K, respectively (see Table 1). This large T c decrease is 
to be ascribed to the strontium ionic radius which is smaller than the barium 
radius. Since the alkaly-earth ions are placed at the centre of the perovskitic 
cells, the smaller Sr +1 ion radius reduces the orbital overlap and thus pairing 
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energy Wp. In Table 2, using the edge, face, centre [E F C] notation [30, 46], 
the probable layering scheme of Bao.6Ko.4Bi03 is shown together with those of 
the previously mentioned BaPbo.7Bio.3O3 and SrTi03_,5 compounds. The ion 
arrangement which originates superconductivity is the same in all compounds. 
In particular, in all compounds an empty C-position separates the alkaly-earth 
ions, thus allowing formation of bonds as occurs in bilayered cuprates. 



[O - Sr] 
[Ti0 2 -] 

f [ Sr+1 ] } 

[Ti0 2 -] 

I [ Sr* 1 ] J 

[Ti0 2 -] 
[O - Sr] 
SrTi0 3 - 5 
T c = 0.3 K 



[O - Ba] 
[Pb0 2 -] 

( [0-BB+ 1 ] ) 

{ [Bi0 2 -] 

I [0-BB+ 1 ] J 
[Bi0 2 -] 
[O - Ba] 
[Pb0 2 -] 
BaPbo.7Bio.3O3 
T c = 10 K 



[Bi0 2 -] 

[O-K] 

[Bi0 2 -] 
f [O-Ba^] I 
4 [Bi0 2 -] 
I [O-Ba* 1 ] J 

[Bi0 2 -] 

[O-K] 

[Bi0 2 -] 

[O-K] 
Bao. 6 Ko.4Bi0 3 
T c = 30 K 





[-0 2 Cu] 




[Ba - 0] 


[-0 2 Cu] 


[0 - Tl] 


[La - 0] 


[Tl - 0] 


f [0-Sr+ 3 ] 1 


f [--Ba+i] 


I [Cu0 2 -] J 


1 [Cu0 2 -] 


[O - La] 


[0 - Ba] 


[La - O] 


[Tl - 0] 


[-0 2 Cu] 


[0 - Tl] 


Lai.85Sro.i 5 Cu0 4 


[Ba - 0] 


T c = 35 K 


[-0 2 Cu] 




Tl 2 Ba 2 Cu0 6 _ 5 




T c = 85 K 



Table 2: Probable layering schemes in mixed stoichiometry perovskitcs. 
Braces enclose the layers which activate superconductivity. Ions with unpaired 
electrons are represented by bold type. 

Probable layering schemes of Lai.85Sro.isCu04 and Tl 2 Ba 2 Cu06-5 com- 
pounds are also shown in Table 2. These cuprates are of special interest since 
both are characterized by isolated Cu0 2 layers. In the oxidized Lai. 85 Sr .i5CuO4 
compound, trivalent Sr +3 ions are included showing unpaired electrons in the 
krypton shell and forming bonds with copper. The situation, however, is quite 
different from that of bilayered cuprates. In fact, while copper is placed in 
E-position, strontium lies in the C-position of the overhanging layer. Conse- 
quently, each copper is allowed to form bonds with four strontium ions so that 
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the symmetry between the interacting layers peculiar to bilayered cuprates no 
longer exists. A possible equivalent interpretation assumes that divalent stron- 
tium causes an equal number of lanthanum ions to be oxidized to valence four, 
thus showing unpaired electrons in the xenon shell. The monolayered reduced 
compound Ti2Ba2Cu06-<5 shows a similar situation. In this compound, mono- 
valent Ba +1 ions with unpaired 6s electrons are present forming bonds with 
3d-electrons of divalent copper. Both the compounds dealth with are indeed 
characterized by staggered overlaps of the interacting layers ( ). This would 
require some modifications to calculations of Section 4, leaving however the es- 
sential results unchanged. In reality, the matters presented in this Section are 
based in part on conjectures owing to the lack of data on the actual placing of 
the unpaired electrons. Notwithstanding this, in our opinion, the reliability of 
the interacting layer mechanism is reasonably proved. 

8. Discussion and conclusions 

According to the electron pairing mechanism we propose, superconductivity 
in cuprates requires the presence of two neighbouring CuC-2 layers. In com- 
pounds like YBa2Cu3C>7, each couple of layers constitutes an independent super- 
conductor. Contrary to this point of view, evidence has been claimed for nonex- 
istence of superconductivity in an isolated Cu0 2 bilayer. Organic chains (Py- 
C n H 2n +i)2Hgl4 (2<n<12) were intercalated in the bilayered Bi 2 Sr 2 CaCu208 
compound thus drastically increasing the distance between consecutive bilayers 
[47]. In this way, a complete disappearance of superconductivity was observed. 
This result was considered as a proof that superconductivity depends on a three- 
dimensional linkage between the couples of neighbouring CuC>2 layers. In oppo- 
sition to this conclusion, we point out that pyridine is an efficient electron donor 
(see for instance [48]). Consequently, the observed disappearance of supercon- 
ductivity is an expected donor-effect similar to that of yttrium-praseodimium 
substitution in YBa2Cu30y discussed in item 3) of Section 6. 

The question of the actual number of layers required for originating su- 
perconductivity in cuprates is certainly of primary importance. The idea that 
interlayer coupling plays a role in superconductivity dates from 1987 when Z. 
Tcsanovic proposed a mechanism that involves Coulomb interaction from the 
band at the Fermi surface to some of the fully occupied or empty bands away 
from the Fermi level [49]. Another momentous question is the basic interaction 
which allows the formation of electron pairs. In alternative to the phonon cou- 
pling peculiar to the BCS theory, in 1997 we proposed the inter-layer HL-type 
two-electron exchange [2] , while T.M. Mishonov ct al. proposed an inter-atomic 
two-electron exchange [50]. In contrast with this previous proposals, these au- 
thors recently have advanced the intra-atomic exchange of two electrons between 
4s and 3& x 2_ y 2 states as the origin of high T c superconduction in cuprates [51]. 
In our opinion, identification of exchange interactions as the very cause of su- 
perconductivity represents a major progress in this field of studies. But, for a 
full understanding of the phenomenon, the question remains to be settled if an 
unique mechanism or different mechanisms are active in the different kinds of 

11 ) A state of affairs of this kind has been already considered in Ref. [3]. 
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unconventional superconductors . 

Leaving aside the theoretical issues, the most pressing thing appears to be the 
discovery of new superconductors of higher T c and, hopefully, of a 300 K super- 
conductor. In reality, even a minor increase of pairing energy Wp may originate 
a large increase of T c , owing to the exponential dependence of Ao on Wp given in 
Eq. . A sound argument in favour of this expectation is offered by the detec- 
tion of a sharp superconductive transition at 235 K due to traces of an unidenti- 
fyed phase fortuitously mixed to a HgBa 2 Ca2Cu30s sample [52] . The approach 
we advise for improving critical temperature is based on the fact that all cuprates 
with divalent calcium sandwiched between the neighbouring CuC>2 layers show 
values of T c higher than 92 K, the YBa2Cu30y critical temperature, in which 
trivalent yttrium is sandwiched between the layers. This fact can be explained 
by considering that the charge of unpaired electrons on copper ions is a little 
shifted towards the sandwiched positive ions, thus reducing overlap of unpaired 
electrons and, consequently, Wp. Obviously, this effect is expected to be less 
harmful with divalent calcium than with trivalent yttrium. This induces us to 
consider compounds of stoichiometry M +1 M^" 3 Cu307 derived from YBa2Cu30y 
by substituting Y +3 with monovalent M +1 ions and Ba +2 with trivalent M +3 
ions (M +1 = Li, Na, K; M +3 = Y, La). This substitution leaves the cell neutral- 
ity unchanged. Using the [E, F, C] notation [8, 9], the substitution spoken of 
is shown in Table 3. Along the same line of reasoning, in the HgBa2Ca 2 Cu30g 
compound we can consider the substitution of Ca +2 ions with M +1 ions and 
Ba +2 ions with M+ 3 ions yielding the HgM^lVLj" 1 Cu 3 8 compound. Also the 
thallium based compounds Tl2Ba2CaCu20s and Tl2Ba2Ca2Cu30io are in prin- 
ciple suitable for substitution of calcium with monovalent ions. Of course, the 
possibility of obtaining these substituted compounds is a mere conjecture which 
should be confirmed by experiments. 



[CuO -] 
[O - Ba] 



[CuO -] 
[O - M+ 3 ] 




[O - Ba] 
[CuO -] 



[O - M+ 3 ] 
[CuO -] 



YBa 2 Cu 3 07 
T c = 92 K 



M +1 M+ 3 Cu 3 7 
T c =?K 
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[Hg- -] 


[0 - Ba] 


' [Cu0 2 -] ' 




[- - Ca] 




< [Cu0 2 -] 




[- - Ca] 




k [Cu0 2 -] , 




[0 - Ba] 


[Hg - -] 



T+31 



[Hg 

[O - M" 1 

[Cu0 2 -] J 
[- -M+ 1 ] 

[Cu0 2 -] 
[- -M+ 1 ] 

[Cu0 2 -] 
[O - M+ 3 ] 

[Hg - -] 



HgBa 2 Ca 2 Cu 3 08 
T c = 133 K 



HgM^M^CusOg 
T, =? K 



Table 3: Layering schemes of YBa 2 Cu 3 07 and HgBa 2 Ca2Cu308 cuprates and 
their respective modified counterparts M +1 M 2 f3 Cu 3 07 and HgM 2 f3 M 2 fl Cu 3 8 . 
Braces enclose the layers which activate superconductivity. 



Appendix A - Some remarks about the Belyaev Hamil- 
tonian for pairing correlations in fermion systems 

S.T. Belyaev during a stay at the Institute for Theoretical Physics of the 
University of Copenhagen wrote a paper entitled " Effect of pairing correlations 
on nuclear properties". We report hereinafter Subsection 1 (pag. 7) of this 
paper ( 12 ) in which the Hamiltonian for the system of interacting particles is 
given. 

oooooooooooooooooooooooooooo 

1. Hamiltonian 

We consider a system of nucleons which are moving in a certain axially 
symmetric self-consistent well. (For simplicity we do not distinguish between 
neutrons and protons). As basic functions of the second quantization represen- 
tation we choose the wave functions of a nucleon in this well. States which differ 
only in the sign of the projections of angular momentum along the symmetry 
axis are degenerate. We call such states "conjugate" states and mark them with 
the index ka = (k+; fc— )* . 

The wave functions of the conjugate states are assumed to transform into 
each other by complex conjugation and exchange of the spinor components**. 

12 ) Published in Matcmatisk-fysiskc Mcddclelser (31, no. 11, 1959) a journal issued by the 
Royal Danish Academy of Sciences and Letters. 
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Let us introduce the Fermi operators a ka ; a ka which create and destroy a 
particle in the state ka. The Hamiltonian for the system of interacting particles 
is then 

H ' = J2 £ k i a t+ a k+ + at_a k _) 

k 

-\ ( k i a ^ a z\ G \ k W2K°\)4 1 «A2V2 a k> 2 *> 2 a K<ji 1 i (!) 

(fe, a) 

where is the single-particle energy in fc-th state. (The sign of G is chosen to be 
positive for an attractive interaction). The Hamiltonian (1) describes a system 
with a fixed number of particles N . Therefore, in a perturbation treatment in 
which H' is split into two parts, each of these parts must commute with N. 
The problem is essentially simplified if we make a transition from the system 
with fixed N (" iV-system" ) to one with a fixed value of the chemical potential 
A (" A-system" ) , which is described by the Hamiltonian 

H = H' — XN. (2) 

The choice of A determines only the average value of N in the A-system. There- 
fore, the solution which corresponds to the Hamiltonian, (2), will describe only 
average properties of nuclei and does not pretend to describe the individual nu- 
clear properties for which one needs a fixed value of N. As will be shown later, 
the uncertainty in the value of N is small. In practice, the averaging is done 
only over a few neighbouring nuclei, either all even or all odd. 

Footnotes 

*) In fact, even symmetry of reflection in a plane is enough for the 
definition of the conjugated states. We speak of axial symmetry only 
for definitencss. 

**) If = (fy), then V>_ = The transformation ^+ -» V>- 

is equivalent to the time reversal T. 

oooooooooooooooooooooooooooo 

Let the Hamiltonian for pair formation, given in Eq.(l) of the previous 
Belyaev's paper, be written in the form 

(fe, a) 

By appropriately defining indexes k and a, it assumes different meanings. In- 
deed, by putting 

fci = k, k 2 = —k, k[ = k', k' 2 = —k', o\ = a\ =|, a 2 = <j' 2 =1, 

(39) 
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\ {k 1 a x k 2 n 2 \G\k'^ l 2 k' x o\) = ^k\-k\\G \-k' { k' |) - V*, fc (40) 
and moving toward left the apex (/), H pair can be rewritten in the quick form 

H P air = - V k, k' Sfc/ T S+ fe , i S_ fei a feT , (41) 
k, k' 

where arrows stand for spin components. This is the Hamiltonian for the Cooper 
pairs in the BCS theory. These pairs are characterized by opposite electron 
momenta, that is, by reflection symmetry around a center. In a different way, 
by allowing for the conjugated states considered in the Belyaev's paper, we have 

fci = k 2 = k, k[ = k' 2 = k' <7i = a\ = +, a 2 = (r' 2 = -, (42) 

^(k 1 a 1 k 2 a 2 \G\k 2 a' 2 k' 1 <T' 1 ) = X - (k+, k-\ G \k'-, k'+) = V k ,, k (43) 

and 

H P air = - X] v k, k> a£, + a£,_a k _a k+ (44) 

fe, k' 

which is the Hamiltonian for the Majorana pairs, characterized by opposite signs 
of the projections of angular momenta along the symmetry axis. We emphasize, 
finally, that as pointed out in footnote (*) of the Belyaev's paper, even fermions 
of equal momenta which exibit symmetry of reflection in a plane can be paired 
in conjugated states. The corresponding Hamiltonian can be written as 

Hp air = - ^2 V k , k' al, R al, L a kL a kR , (45) 

k, k 1 

where indexes L and R mean left and right with respect to the reflection plane. 
Symmetries of interacting fermion pairs are summarized in Table 4. Pairs show- 
ing planar symmetry have been named " Lewis pairs" . This because, as explained 
in Sections 4 and 5, superconductivity in cuprates can be ascribed to pairs of 
electrons forming covalent bonds characterized just by this kind of symmetry. 

Table 4: Symmetries of fermion pairs. 

Pairs Conjugated states Symmetry 

Cooper —k,l * + k, | Central 
Majorana k, — * k,+ Axial 
Lewis k, L * k,R Planar 
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Appendix B - The electron energy spectrum 

As shown in Eq. © and by writing briefly Wk for W a (k ) and Wb (k D ), 
the energy of electrons running on the C11O2 layers is given by the expectation 
value 



1 / N 

Wk = 77 ( X! cxp k ' u ™) a ( r ~~ Um ) 



of Hamiltonian 



H(p,r) 



^2 ex P (i k • u„) a (r - u n )j 

(46) 



n=l 



A? 



Ze 2 



+ ^a* (r) . 



(47) 



0=1 



This energy is of course the same for a or b layers. For completeness sake, in 
Eq. (|47|) the lattice potential Vi at (r) has been included due to ions in positions 
other than u p , that is, to oxygen, yttrium and barium ions neighbouring the 
copper ions. Orbitals of copper ions are solutions of equation 



2m 



Ze 2 



+ Viat (r) 



a (r - u) = E 3d a (r - u) 



(48) 



Ezd standing for the orbital energy. On the CuC>2 layers, copper ions form a 
square grid with Cu +2 ions at the square vertices and 0~ 2 ions at the middle of 
the square sides. The Coulomb field of O -2 ions cuts down the electron charge 
density of Cu +2 ions along the square sides, thus increasing density along the 
square diagonals. So the electron charge distribution is expected to show the 
four- lobe shape peculiar to d-orbitals ( 13 ). 

Cu +2 ions placed at opposite ends of the square sides cannot interact directly 
owing to the interposite oxygens. On the contrary, Cu +2 ions placed at opposite 
ends of the square diagonals are able to originate overlap and exchange integrals 
of significant values. This entails that Cu +2 ions placed alternatively along the 
square sides form two independent but equivalent ion sets, each holding N/2 
ions. It is therefore sufficient to consider one of these sets. So Eq. can be 
rewritten as 



W k = — ^2 exp [i k • (u„ - u m )] / a (r - u m ) H (p, r) a (r - u„) d 3 r. 

m, n— 1 

(49) 

By taking into account that for n 7^ m only ions in the four neighbouring lattice 
positions Ui around u m make a significant contribution, we have 

13 ) As for the real form of copper orbitals in C11O2 layers, the most likely assumption is that 
they consist of a superposition of 3d xy and 3^2 _ y 2 orbitals. The 3d x 2 _ y i orbitals are lined 
up along the Cu-O-Cu chains, thus allowing for super-exchange interactions between coppers 
mediated by oxygen ions. Consequently, only the ?>d xy contributions should be considered 
when dealing with pairing energy of electrons on neighbouring Cu02 layers (see Ref . [4] ). 
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N 



N/2 r 

E 

rn — 1 



a (r m ) (p, r m ) a (r m ) d 3 r m + 



^ exp (£k • <Xi) / a (fj + <x;) i? (p, r,) a (r<) d 3 r; 



(50) 



where r TO = r u m , = r and <Xj = u m . For an unlimited CuC>2 
layer, the sum over m is independent of position u m so that all terms in the 
sum are equal. This allows us to write 

W k = Ja{r) H(p,r) a (r) d 3 r+ 
+ ^exp (ik • o-j) J o(r 4 + <t;) F (p, rj a(r i )d 3 r i . (51) 



On the other hand, we have from Eqs. I|48|) and l|47|l 



Ze 2 



H (p, r m ) a (r m ) = I E 3d - ^ — a (r m ) , 



(52) 



in which, as for Eq. (|50|l . only terms for ions in the four positions around 
u m have been included. Owing to square symmetry of these positions around 
u m , we obtain 



where 



J a (r) H (p, r) a (r) d 3 r = E 3d - 4£ c , 

E c = a r 1 , a r d 3 r 

7 |r-cr| 



(53) 



(54) 



means Coulomb energy of the ion in position u originated by the electric field 
of the ion in position u + cr. By putting r p = r — u p , we have analogously 

y a (r, + <7j) if (p, rj a (r*) d 3 r i= 

/4 /. Z e 2 

a(r, ; +o-j) a(ri)d 3 ri - ^ / a ( r i + °») "1 — r a ( r <) d3r .- ( 55 ) 

Since the product a (r^ + <Tj) a (r,) shows a significant value only midway be- 
tween the u m and Ui positions, only the integral for p — m is to be accounted 
for. By omitting label i and introducing the overlap 



= J a (r + cr) a (r) d 3 r 



(56) 
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and exchange 



E ex = fa(r + a) a (r) d 3 r (57) 

integrals, we obtain 

J a ( Ti + m) H (p, r t ) a (r 2 ) d 5 v t =OE M - E ex . (58) 
In this way, substitution of Eqs. I|53|) and l|58[l into Eq. 151|) leads to 

4 

W fe = E 3d - AE C + (OE 3d -E ex )J2 exp (ik • O . (59) 

t=i 

This result mimics the one obtained in the case of tridimensional lattices in 
which six vectors (Ti have to be considered. 

By assuming x and y axes parallel to vectors (Ti, by letting Wo = E 3 d—^Ec J r 
4 (OE 3 d ~ E ex ) and B = 4 (OE 3 d — E ex ) for ground state and band energies, 
respectively, the electron kinetic energy turns out to be 

£fe = Wfc - Wo = — [2 - cos (fc x cr) - cos (fcj,er)] , (60) 

where a = |<Tj| . It follows that for fc x er, k y a = ±7r, kinetic energy attains its 
maximum value 2B. By assuming the CuC>2 plane to be a square of area a x a 
with sides parallel to x and y axes and taking into account that k x = iij (tv/cl) , 
fc. y = n y (jr/a) with n x , n y = 0, ±1, ±2..., Eq. H6()|l can be rewritten as 



n x n y \ B 



n n J 2 



n x \ / n y 

7T — COS IT — 

n J \ n Q 



(61) 



where no = a/ a. To find the isoenergetic contours on the n x , n y plane, the 
initial values n x /no = and n y /no — are chosen corresponding to energy 

e(0) = f [1-cos(tt0)]. (62) 

For 6 = 1, we have e(l) = B which is the maximum value of kinetic energy 
on the isoenergetic contours. This means that B represents the actual band- 
width. Then, by keeping e(n x /no, n y /no) = n y /no is evaluated as a 
function of n x /no for various values of 9 by means of a numerical procedure 
which also finds the area enclosed in the contours. Owing to electron spin, 
twice this area represents the number N s of states of energy less than e (9) . 
For 9=1, the contour is a square of half-diagonal no and area 2tt,q (see Figure 
1). Thus, for 9 = 1, N s assumes its maximum value N s q = 4riQ. On the other 
hand on the Cu02 planes each mesh of area cr 2 /2 holds one electron, so that the 
overall number of electrons is 2uq. This means that the band is half-filled. For 
N s /N s o = 0.5, we find 9 = 0.59. It follows, from Eq. (JHU, s F = 0MB or, by 
identifying the Fermi energy with the chemical potential, a = 0.64£>. Utilizing 
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the previously mentioned numerical data, the density of states per unit cell at 
the Fermi level is found to be 

n F = l^£ = l^, (63) 
N deF N B y ' 

which, taking into account that 2n§ = N, leads to the simple relationship 

/j,fl F = 1.28. (64) 

It is to be pointed out that this result holds in general independently of band- 
width B, that is, of the actual values of overlap and exchange integrals. 



Figure 1: Isoenergetic contours for electrons running on CuC>2 planes. The 
square contour for n x /n = 1 corresponds to band- width energy B, the one for 
n x /n — 0.59 to the Fermi energy ep- 



Appendix C - 

systems 



Conservation of momentum in fermion 
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Hamiltonian i|21[l conserves momenta of the interacting electrons. To realize 
how this occurs, we recall that in degenerate fermion systems, by letting n r 
be the number of particles on the r-th level and u) r the level degeneracy, we 
have, at zero temperature, n r = tu r below Fermi level and n r — above it. 
Consequently, taking into account that fluctuation of n r is 

5n r = \/n r {l — n r /uj r ), (65) 

no fluctuation is allowed [53]. For the same reason, no scattering is possible 
either ( 14 ). This is a direct consequence of the Pauli principle. A clear evidence 
about this property of degenerate fermion systems is provided by the nuclear 
matter. In fact, nuclei can be described both by the liquid drop model, in which 
nucleons strongly interact among themselves, and by the shell model, in which 
nucleons behave as particles moving freely in a potential well. This confirms 
that in degenerate fermion systems interacting particles maintain unchanged 
their momenta. 

We must therefore verify that Hamiltonian (|21|l really does not allow for 
scattering processes. We have 

Hp a i r = -W P ^ ^k'a^k'b^kb^ka- (66) 
k, k' 

Terms with k ^ k' change pairs with momenta k, a and k, b in pairs with mo- 
menta k',a and k',b. But each term with k — p and k' = q is associated 
with a term with k = q and k' — p. These terms represent opposite scattering 
processes. Omitting factor —Wp, their contribution in H pa i r is 

Spq = Uqa&qb&pbapa + ap a a pb a qb a qa . (67) 

Taking into account the Jordan and Wigner anticommutation rules, operator 
Sp q can be rewritten in the ordered form 

Spq = a+ a a pa a+ b a pb + a qa a pa a qb a pb . (68) 

Operators a and a + and fermion states |0) and |1) can be represented by the 
matrices 






1 







, |o> = 


1 


, ll> = 













1 





1 



(69) 



We have, indeed, 

a|l} = |0), a + |0) = |l), a|0) = 0, a + 11) = 0. (70) 

14 ) For T > K, states in the Fermi level region are only partially occupied. Thus, in this 
region scattering processes can take place. But it is wrong to ascribe superconductivity to 
these scatterings because when T goes to zero scatterings are shut off while the supercon- 
ductig energy gap Ao attains its maximum amplitude. Fermi level appears in Eq. 1321 since 
A is small with respect to chemical potential fj,, as explained in Eq. 1301 . 
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We consider, moreover, the hcrmitian operator 



X = a + a = 
which interchanges the fermion states, that is 



1 

1 



x|o) 



1 



2|i) = |o) 



(71) 



(72) 



In this way, operator S pq can be written as the sum of two 8x8 matrices, each 
showing on the diagonal four 2x2 matrices for the a and a + operators cor- 
responding to qa, pa, qb, pb states. So, by taking into account Eq. 1|71JI . we 
have 



^PQ ^Cqa\pa^(-qb\pb. 

By considering the four fermion state |1, 1, 1, 1) = 
tain 



\i) qa m pa w qb w pb 



;i,i,i,i|s* P9 |i, 1,1,1) = <i,i,i,i| o,o,o,o) = o, 



(73) 
we ob- 

(74) 



that is, the expectation value of S pq vanishes. The same result is found in general 
for all fermion states. Consequently, no scattering process is really accounted 
for in Hpair Hamiltonian. Only non-scattering terms with k = k! give a contri- 
bution. This means that electron momenta are conserved in Lewis pairs. This 
argument applies unchanged to Cooper and Majorana Hamiltonians in which 
scattering terms can be eliminated in the same way. The lack of effective scatter- 
ing processes entails that momentum of fermions is conserved for all the pairing 
mechanisms considered in Table 4. 
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